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Аннотация 

Введение. Основу для анализа и интерпретации данных гидродинамических исследований составляют 
теоретические модели и вычислительные алгоритмы. Но, несмотря на востребованность данной тематики, 
многие вопросы, связанные с нестационарными течениями жидкостей в нефтяных пластах, все еще требуют 
решений. Поэтому математические постановки задач, связанных с учетом нестационарного течения жидкости, 
разработка эффективных численных методов и алгоритмов, их решение с применением современных 
уеБ-технологий являются актуальными. Цель данной работы — разработка \еб-приложения для 
математического моделирования процесса фильтрации жидкости в однопоровых коллекторах при проведении 
гидродинамического исследования на добывающей скважине. 

Материалы и методы. Для решения поставленной задачи применены методы механики сплошной среды 
и вычислительной математики. Представлена модель течения нефти в однопоровом коллекторе. При разработке 
приложения использованы языки программирования Руфоп и Лауабсире. Результаты расчетов хранятся 
в реляционной базе данных, реализованной средствами РозёотеЗ ОГ. 

Результаты исследования. Разработано новое \еБ-приложение для моделирования процесса фильтрации 
нефти в однопоровых коллекторах, которое применимо для изучения флюидодинамических процессов и может 
быть использовано для прогноза дебитов, добычи и расчета оптимальных режимов работы скважин. 
Обсуждение и заключение. Разработанное у\еб-приложение позволяет строить поля давления и температуры 
в пласте около работающей и остановленной добывающей скважины и на различных расстояниях от нее. 
Данная информация дает возможность оперативно оценивать длительность гидродинамических исследований, 
а также регулировать работу скважин. Приложение может быть развернуто в существующей сетевой 
инфраструктуре, пользоваться всем функционалом, подключаясь к удаленному серверу. Оно оптимизировано 
для использования на различных платформах и имеет широкие перспективы дальнейшего развития. 


Ключевые слова: у\еБ-приложение, математическая модель, гидродинамические исследования, база данных, 
давление, температура 
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Введение. Российская Федерация является одной из ведущих стран мира по объемам производства нефти. 
Каждый год на месторождениях страны добывается свыше 500 миллионов тонн нефти, подавляющая часть 
которой в сыром виде экспортируется в страны азиатско-тихоокеанского региона, Европы и ближнего 
зарубежья. Развитие нефтяной промышленности в нашей стране оказывает влияние на другие отрасли и сферы 
деятельности, поэтому решение задач по эффективному освоению и разработке нефтяных месторождений было 
и остается перспективным. Основой современных научных технологий нефтяной разработки является 
комплексное изучение свойств продуктивных коллекторов, содержащихся в них веществ и происходящих 
сложных процессов при течении флюида. 

Главная цель при добыче нефти — увеличение темпа отбора флюида из продуктивного пласта 
и поддержание пластового давления. Планирование работ по эффективному освоению и оптимальной 
эксплуатации нефтяной скважины носит комплексный характер, связанный с гидродинамическими процессами, 
происходящими при течении жидкости в системе «скважина — пласт» [1]. Для изучения продуктивных свойств 
пластов и жидкостей применяются гидродинамические исследования скважин (ГДИС). Средствами ГДИС 
проводятся замер совокупности данных коллектора с использованием глубинных или устьевых приборов, 
атакже последующая обработка замеров и интерпретация полученных данных (продуктивность, 
фильтрационные свойства жидкостей, тип коллектора и т. д.) [2]. 

При исследовании коллекторов с использованием ГДИС применяются теоретические модели, основанные 
на классических уравнениях и законах, а также вычислительные алгоритмы для моделирования процессов 
фильтрации флюидов. Однако, несмотря на востребованность данной тематики, многие вопросы, связанные 
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с математическими и вычислительными аспектами моделирования нестационарных потоков в коллекторах, 
требуют дальнейшей разработки и развития [3, 4]. 

Для анализа и интерпретации результатов ГДИС существует ряд отечественных («Автон», «Гидрозонд») 
и зарубежных (ОСА, ГедаРю\м, Но\мУ1910оп, МАЗТ, барыг МГ.) программ [5-7]. При всех имеющихся 
преимуществах они обладают некоторыми недостатками, такими как невозможность проведения полного 
спектра расчетов, просмотра и изменения параметров исходных моделей. Все представленные программные 
продукты являются настольными, коммерческими, сложными в эксплуатации, многие обладают 
перегруженным интерфейсом. В условиях, когда необходимо оперативно или в режиме реального времени 
провести расчет, большие коммерческие симуляторы оказываются малоэффективными. Поэтому для 
выполнения схожих работ предлагается использовать \еБ-приложение, обладающее множеством преимуществ. 

В отличие от настольных, У\еЬ-приложения независимы от операционной системы и вычислительных 
мощностей персональных компьютеров. Программный код \еЬ-приложения пишется единожды для 
конкретной платформы, где оно будет развернуто. У!еБ-приложения выделяются простой реализацией, 
легкостью применения, высокой производительностью. Их основными характеристиками являются 
доступность, надежность, безопасность, масштабируемость, гибкость, кросс-платформенность. Для работы 
с \еЬ-приложением от персонального компьютера пользователя необходимо находиться в одной сети 
с сервером, на котором оно развернуто. Такой сетью может быть локальная сеть предприятия или глобальная 
сеть Интернет, если речь идет о веб-сайтах. По сравнению с настольными, у \еБ-приложений можно выделить 
один недостаток — время отклика, зависящее от удаленности персонального компьютера от сервера. 

Логика работы \еб-приложения отличается от настольного, которое запускается и выполняет код 
на компьютере пользователя. Поэтому настольные приложения отличаются более расширенным и отзывчивым 
пользовательским интерфейсом и позволяют реализовывать более сложную бизнес-логику работы. Однако 
с развитием \уеб-технологий это достоинство с каждым годом нивелируется, и \еб-приложения все чаще 
бывают способны реализовывать сложные пользовательские сценарии. Таким образом, в части реализации 
производственных задач \еб-приложения становятся более конкурентными. 

В связи с этим математические постановки задач с учетом нестационарного течения жидкости во всех его 
элементах (пласт, стволовые трубы, проточные каналы), разработка эффективных численных методов 
и алгоритмов, последующие их решения с применением современных технологий \е-разработки являются 
актуальными. Использование полученных численных алгоритмов позволяет количественно описывать 
поведение давления и температуры в нефтяных скважинах. Алгоритм и \еБ-приложение предоставят 
пользователям, специалистам в области анализа и интерпретации гидродинамических исследований, удобный 
и легко адаптируемый к их практическим потребностям инструмент для расчета сложных процессов, 
происходящих при разработке пористых коллекторов. 

Материалы и методы. Рассматривается добывающая скважина, на которой проводится гидродинамическое 
исследование методом кривой восстановления давления. Скважина эксплуатирует терригенный коллектор, 
который относится к однопоровым. Пласт является неограниченным, воздействие соседних скважин 
отсутствует. Течение флюида в пласте — горизонтальное, перетоки между пропластками отсутствуют. 
Скважина работает в установившемся режиме с постоянным забойным давлением. На границе контура 
поддерживаются постоянные давление и температура. Для проведения исследования измерительный прибор 
опускается на глубину верхних дыр перфорации. Расчет динамики давления и температуры проводится 
в пространственно-одномерном случае. 

При численном моделировании важно учитывать все параметры пласта и флюида, чтобы достичь наиболее 
качественного описания процесса массопереноса нефти в однопоровом коллекторе. Рассмотрим 
математическую модель, описывающую процесс фильтрации нефти в однопоровом коллекторе. Основой для 
модели служат уравнения пьезопроводности и теплопроводности. Математическая модель в радиальных 
координатах представляет собой систему дифференциальных уравнений в частных производных: 
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гдеф — пористость пласта; с, — сжимаемость породы (1/Па); к — проницаемость пласта (м); и — вязкость 


жидкости (Па с); Р — пластовое давление (Па); г — расстояние до стенки скважины (м); # — время (с); с — 
теплоемкость нефти (Дж/кг К); р — плотность нефти (кг/мЗ); Т — пластовая температура (К); А — 
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теплопроводность пористой среды (Вт/м?-К); и — скорость конвективного переноса тепла в пористой 
среде (м/с); = — коэффициент Джоуля-Томсона (КЛЛа); п — коэффициент  адиабатического 
расширения (К/Па). 
Первое уравнение описывает процесс изменения давления в пласте, второе — процесс теплопереноса жидкости. 
Скорость конвективного переноса тепла в пористой среде определяется выражением 


КОР 
и=-——. (2) 
и д" 
Система уравнений (1) задается в пространственно-временном интервале: 
Г, <Г< К, (3) 
О<1<4, 
где г, — радиус скважины; № — расстояние до границы контура исследования (м); #, — время исследования (с). 


Исходя из поставленных выше условий для скважины и рассматриваемой области для системы (1) 
начальные и граничные условия задаются следующим образом: 
1=0:Р=А,Т =Т,г, <г< К, 


г=г,:РЕВ,Г=Т,1>0, (4) 
т и о 0,Е> 0. 
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После определения начальных и граничных условий получена модель процесса фильтрации нефти в 
однопоровом коллекторе (1)—(4), характеризующая перераспределение давления и температуры в пласте. 

Первоначальным этапом решения системы является приведение его к дискретному виду. Для этого 
используется метод конечных разностей [8]. Строится равномерная пространственно-временная сетка: 
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(5) 
= 1 
—_ И —_ — к 
С, = |. =т,п = 0,1,..., М, = т 
где й — Шаг сетки по пространству; т — шаг сетки по времени; № —_ число узлов по пространству; М — 
число узлов по времени. 
После замены дифференциальных выражений разностными аналогами получается система уравнений: 
р"и т Чт р"и Чт р" 
р" Р" | 1 Е р па п ? ги — 
фс о 2 2 5 
ф ) 
< ых п? 
п+1 п п+1 п+1 п+1 п+1 п п+1 п+1 
ТТИ АИТ ОТ Ти пы ВЫ". ь 
ср : сри, 8 (6) 
т Г р т й 
р"и _ рп 
фер 1 =2,...,М-1п> 0. 
т 


Система (6) представляет собой систему линейных алгебраических уравнений (СЛАУ), которая сводится к 
трехточечному уравнению общего вида: 


Аруги — Ву" + Су =ЕМ=1,...М,п>0. (7) 
Коэффициенты равны: 
т ТГ Г 
АР - — ‚В; — - 2 т ме — — в = ыы В) 
НИ т НГ т НА п т 
^, 2^ сри;, ср ^, 
А’ 2.0 ‚Ву =| 72 + ме. 2.’ (8) 
Й`т йух т т Й-т, 
п+1 п+1 п+1 п 
”" -Р Р"` -Р 
ЕТ = "(и +1) 7 +срив- > : 1 фср : 


Система (7) решается методом скалярной прогонки [9, 10]. Задача является нестационарной, поэтому метод 
применяется для определения пространственного распределения давления и температуры пласта на каждом 
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временном слое [11, 12]. Задача решается последовательно, на первом этапе решается уравнение для давления, 
на втором — для температуры. 

Для автоматизации решения задачи разработано \еб-приложение. В общем виде его структура представлена 
на рис. 1. В ней выделяются пять основных компонентов: входные информационные потоки, выходные 
информационные потоки, методы обработки информации, технические средства обработки информации, база 
данных [13, 14]. 


Методы обработки Технические 
информации средства обработки 
информации 


Входные Выходные 
информационные У’еБ-приложение информационные 
потоки потоки 


База данных 


Рис. 1. Структура приложения в общем виде 


Входными информационными потоками являются данные о параметрах модели, а выходными — графики, 
отражающие временные или пространственные зависимости забойного давления и температуры. Под методами 
обработки информации понимаются численные методы и методы интерпретации гидродинамических 
исследований. Техническими средствами обработки информации являются сервер, на котором размещено 
приложение, и персональный компьютер пользователя. Для хранения входных параметров и результатов 
расчетов используется база данных. 

Конечный вид программного обеспечения зависит от средств, выбранных для его разработки. Средствами 
разработки являются парадигмы и языки программирования, среды и технологии для разработки. Выбор 
инструментария влияет на способ реализации задачи, критерии качества программного обеспечения, 
представление данных и тип графического интерфейса пользователя. Для разработки приложения в качестве 
языка программирования для серверной части выбран Ру#фоп с использованием микрофреймворка Е!азК, для 
клиентской части — ЛауабспрЕ с Уче.]з [15, 16]. В качестве СУБД для проектирования и разработки базы 
данных используется РозотеЗ ОГ. [17]. 

Результаты исследования. Разработанное программное обеспечение представляет собой кросс- 
платформенное клиент-серверное \еБ-приложение. Хранение данных и вычисления происходят на \еЪ- 
сервере, с которым пользователь взаимодействует через браузер (поддерживаются все современные интернет- 
браузеры). Клиентская часть может использоваться на широком спектре платформ (ПК, ноутбуки, смартфоны) 
и операционных систем (\Лпдо\з$, МасО5, Апагола, 105). 

Большую часть главного окна \еб-приложения занимает средство построения графиков температуры и 
давления. Графики строятся после произведения расчета путем выбора соответствующего пункта меню. Также 
имеется возможность построения графиков по ранее произведенным расчетам, для этого выбирается 
необходимый расчет в левой части окна. Графики строятся по пространственным и временным координатам. 
Для изменения вида графика выбираются соответственные пункты в верхней части окна. 

Все произведенные расчеты хранятся в разработанной реляционной базе данных Йо\са|саЬ, структура 
которой представлена на рис. 2. 
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Рис. 2. Структура базы данных Но\’са]саь 


В общей сложности база данных содержит пять таблиц, которые находятся в третьей нормальной 


форме [18]. Все таблицы имеют первичные ключи. Главная таблица, са1сшаНоп, хранит название расчета и дату 


его произведения. Две зависимые от нее таблицы, гаЧга$ и ите, хранят информацию о пространственно- 


временных координатах. Таблицы 


соответственно. 


{етрегаште и ргеззиге хранят значения температуры 


и давления 


В разработанном приложении проведен ряд вычислительных экспериментов, позволяющих оценить 


результат его работы. Для проведения моделирования использованы значения из таблицы 1. 


Таблица 1 


Начальные и граничные условия, параметры расчета 


Параметр Значение Единица измерения 
Радиус скважины, Г, 0,102 м 
Радиус исследования, В 100 м 
Время исследования, 1 24 Ч 
Начальное давление, Ру 25 МПа 
Конечное давление, Рк 30 МПа 
Проницаемость, К 1е-14 м? 
Пористость 0,2 
Вязкость нефти, | 1,10е-3 Па-с 
Сжимаемость нефти 1,58е-9 ИПа 
Теплоемкость, с, 1,4 МДж/ мз-К 
Теплопроводность, А. 0,6 Вт/м-К 
Коэффициент Джоуля-Томсона, = 3,94е-7 К/Па 
Коэффициент адиабатического расширения, 7 3,15е-7 К/Па 
Начальная температура, Т, 363,15 К 
Количество точек по времени, М 1000 
Количество точек по пространству, М 1000 


С использованием данных из таблицы 1 построены графики распределения давления и температуры по 


пространству, представленные соответственно на рис. Зи4. 
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Рис. 3. Распределение давления по пространству в конечный момент времени 
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Рис. 4. Распределение температуры по пространству в конечный момент времени 


Как следует из данных на рис. 3 и 4, после остановки скважины на гидродинамическое исследование 
наблюдается восстановление давления и температуры, причем чем дальше от стенок скважины, тем больше 
величина восстановления. На рис. 5 и 6 показано, как именно распределены давление и температура в 
различные моменты времени (6, 12, 24 ч.) после остановки скважины. Эта информация полезна специалистам 
по анализу и интерпретации гидродинамических исследований, поскольку позволяет оценить динамику 
давления и температуры и принять решение о продолжении или остановке исследования, а также проведении 
дополнительных мероприятий. 
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Рис. 5. Распределение давления по пространству в различные моменты времени 


БИр://уезыиК-доп$.ги 


428 


Мазитов А.А. ИеБ-приложение для математического моделирования нестационарного течения нефти в пористой среде 


Тетрегайхе, К 


——{=24В 
365 ——1=128 


364 —— (=6Ь 


363 


362 
0 20 40 60 80 100 120 


Вад, т 


Рис. 6. Распределение температуры по пространству в различные моменты времени 


На основе произведенных расчетов установлено, что наибольшее влияние на поведение давления 


и температуры оказывает проницаемость пласта. На рис. Тиз изображены графики динамики давления 


и температуры 


в зависимости от различных значений проницаемости вблизи от стенок скважины 


(низкопроницаемый коллектор — 1мД, среднепроницаемый — 10мД, высокопроницаемый — 100 мД 


(1Д = 10`12м?)). Как видно на графиках, при увеличении проницаемости коллектора скорость восстановления 


давления и температуры увеличивается. На рис.7 изображен график для реальных промысловых 
данных (проницаемость коллектора — 80 мД). 
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Рис. 8. Динамика температуры для разных проницаемостей 
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На основе динамики давления можно определить время проведения исследования (рис. 7). Например, для 
высокопроницаемого коллектора не имеет смысла проводить исследование более чем -—6 часов, поскольку 
после этого времени давление практически не меняется, а простой скважины является неэффективным, 
поскольку идут потери добычи. 

Обсуждение и заключение. Разработанное \еб-приложение может быть использовано для анализа 
и интерпретации данных гидродинамических исследований. \еб-приложение обеспечивает удобный 
пользовательский интерфейс и снижает затраты на поддержку. Оно может быть развернуто в существующей 
сетевой инфраструктуре или подключено к удаленному серверу для обеспечения полной функциональности. 
Применение приложения позволяет качественно и количественно оценивать поведение давления и температуры 
в однопоровых коллекторах. А хранение результатов расчетов в базе данных дает возможность проводить 
анализ с использованием результатов множества исследований, на основе которого можно составлять сводный 
отчет и формировать рекомендации по дальнейшему использованию скважины или месторождения в целом. 
Используемые модели являются устойчивыми и сходящимися. Графики строятся с помощью написанного \ер- 
приложения, и поэтому оно может использоваться на производстве. 

Приложение имеет широкие перспективы развития, например, можно реализовать другие модели течения 
жидкости или предусмотреть наличие нескольких фаз. Также на момент написания данной работы остро встал 
вопрос об отказе от зарубежных коммерческих продуктов в пользу отечественных симуляторов. Поэтому 
разработка подобных продуктов является перспективной. 
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